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ABSTRACT 

We present ultraviolet photometry (NUV band, 180-280 nm) of 405 asteroids observed serendipitously by 
the Galaxy Evolution Explorer ( GALEX) from 2003-2012. All asteroids in this sample were detected by 
GALEX at least twice. Unambiguous visible-color-based taxonomic labels (C type versus S type) exist for 315 
of these asteroids; of these, thermal-infrared-based diameters are available for 245. We derive NUV — V color 
using two independent models to predict the visual magnitude V at each NUV-detection epoch. Both V models 
produce NUV —V distributions in which the S types are redder than C types with more than 8er confidence. This 
confirms that the S types’ redder spectral slopes in the visible remain redder than the C types’ into the NUV, this 
redness being consistent with absorption by silica-containing rocks. The GALEX asteroid data confirm earlier 
results from the International Ultraviolet Explorer, which two decades ago produced the only other sizeable set 
of UV asteroid photometry. The GALEX-derived NUV — V data also agree with previously published Hubble 
Space Telescope (HST) UV observations of asteroids 21 Lutetia and 1 Ceres. Both the HST and GALEX data 
indicate that NUV band is less useful than u band for distinguishing subgroups within the greater population 
of visible-color-defined C types (notably, M types and G types). 

Subject headings: surveys — minor planets, asteroids: general — solar system: general 


1. INTRODUCTION 

As in visible wavelengths, ultraviolet flux from asteroids is 
reflected sunlight. However, the steep drop in the solar spec¬ 
trum shortward of ~300 nm (Figure 1) makes asteroids or¬ 
ders of magnitude fainter in the UV than in the visible. For 
this reason—as well as the strong UV absorption by atmo¬ 
spheric ozone—UV observations of asteroids typically em¬ 
ploy the Hubble Space Telescope (HST) or specialized instru¬ 
ments on a space-mission payload physically closer to the 
asteroid. These constraints have generally prohibited large- 
sample demographic studies of asteroids in the UV. 

Predating HST, the International Ultraviolet Explorer 
HUE) targeted 45 asteroids from 1978-1992, producing what 
remains to date the largest published sample of near-UV as¬ 
teroid spectra (Roettger and Buratti 1994), specifically in the 
range of 230-325 nm. The IUE data show evidence of clus¬ 
tering, principally with respect to geometric UV albedo. This 
clustering becomes further evident when coarsely-defined vis¬ 
ible spectral type is included as a categorical parameter for 
each object (C, S and M types being the classes considered 
in the original work). Comparing the /I/E-derived geometric 
UV albedos for each class with the geometric visible albedos 
demonstrated that the S types, which are redder-colored in the 
visible (specifically, 400-800 nm), remain redder than C types 
into the NUV. 

The IUE data, combined with previously-measured visible 
spectra, suggested that asteroid reflectances over the entire 
near-UV to visible wavelength range (200-800 nm) are gen¬ 
erally consistent with those of silica-bearing rocks (e.g., Wag¬ 
ner et al. 1987). To first order, this trend is characterized by an 
increase in a rock’s reflectance at longer wavelengths, gener¬ 
ally attributable to the decreased fraction of volume-scattered 
light, i.e., light which penetrates into the mineral grains. The 
intensity of volume-scattered light varies as exp(— kd/X), 


where d is the grain size, k is the imaginary part of the in¬ 
dex of refraction, and A is the wavelength. Once refracted 
into the grains, volume-scattered light is subject to absorp¬ 
tion by various (A-dependent) interactions with the mineral’s 
crystalline structure. At sufficiently short wavelengths (i.e., 
well into the ultraviolet region), most incident light pene¬ 
trates the grains and is absorbed, while the small amount of 
measured reflected light is predominantly scattered directly 
from the grain surface. Precise characterization of this tran¬ 
sition between surface-dominated and surface plus volume- 
scattered reflectance—as well as the identification of any ad¬ 
ditional mineral-specific spectral features—is therefore useful 
for tying astronomical observations of asteroids to laboratory- 
measured analogs, including lunar and meteoritic samples. 

In this work we aim to verify the IUE’s findings with a 
newer and larger sample of UV asteroid data from the Galaxy 
Evolution Explorer (GALEX), a NASA Small Explorer-class 
space telescope mission which from 2003-2012 conducted a 
UV imaging survey in a far-UV band (FUV, 130-190 nm) 
and a near-UV band (NUV, 180-280 nm). Approximately 2/3 
of the sky was covered, with avoidance of bright stars and 
low galactic latitudes. Martin et al. (2005) discuss the extra- 
galactic science program, while Morissey et al. (2005, 2007) 
discuss the on-orbit performance, survey calibration and data 
products. GALEX has a 50 cm 2 effective area, 1.25 degree 
diameter circular field of view, and full-width-half-maximum 
resolution of 4.5" in the NUV. Programs within the GALEX 
mission included an all-sky survey (AIS, with ~ 100 s expo¬ 
sures) and a medium-depth survey (MIS, with ~ 1500 s ex¬ 
posures), and also a spectroscopic (grism) survey. Figure 1 
shows the photometric response functions of the two GALEX 
bandpasses multiplied by the solar spectrum, with comparison 
to the ugriz visible bandpasses. Detection of asteroids in the 
FUV is extremely unlikely (nonetheless, as described below 
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Table 1 

Observations of asteroids detected in GALEX NUV images. Includes 1,342 detections of 405 asteroids detected at least twice. 


asteroid 

number 

observation 
date (UT) 

detected 
R.A. (deg) 

detected 
Dec. (deg) 

position 
residual ( ;/ ) 

NUV 

mag 

NUV mag 
uncertainty 

MPC-predicted 
visible mag (Vmpc) 

exposure 
time (s) 

unique database ID 
(ob j ID key in CasJobs) 

i 

2011-10-12.65797 

355.22180 

-18.47266 

0.5 

14.38 

0.01 

8.0 

91 

6380556162844065792 

i 

2011-10-21.42095 

353.95617 

-18.39556 

0.2 

14.67 

0.01 

8.1 

80 

6380556163951362048 

3 

2005-12-26.72536 

75.12144 

-1.31299 

1.4 

14.43 

0.01 

7.8 

80 

6381858059773280256 

3 

2011-04-17.04995 

166.80140 

8.52138 

1.4 

16.30 

0.01 

9.8 

1513 

3855329770719936512 

6 

2006-08-29.56208 

309.25038 

-19.31464 

0.9 

14.78 

0.01 

8.4 

112 

6379782093773209600 

6 

2005-05-07.46079 

204.16985 

11.58510 

0.4 

16.67 

0.02 

10.2 

112 

6378656257217134592 

8 

2004-12-21.41462 

122.52225 

19.17684 

0.7 

15.63 

0.02 

9.0 

92 

6377776615736213504 

8 

2004-12-21.48313 

122.51009 

19.18345 

0.1 

15.60 

0.02 

9.0 

87 

6377776615769767936 


This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. 
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Figure 1. GALEX \JX and SDSS visible filter response curves (colored lines) 
multiplied by the spectrum of a G2 V type star (black line). The spectrum is 
from the library of Pickles (1998). The vertical scale is in AB magnitude 
units per unit wavelength, offset by an arbitrary constant. Note wavelength is 
plotted on a log scale. 

we searched for both NUV and FUV asteroid detections). 

Our approach in analyzing GALEX asteroid observations 
differs somewhat from Roettger and Buratti’s treatment of the 
IUE data. First, instead of referencing taxonomic class labels 
( e.g., ‘C type’,‘S type’,‘M type’, etc.) assigned to individual 
asteroids by previous authors, we define classes using a color 
index derived from a clustering analysis performed on a com¬ 
pilation of seven visible-color surveys (Waszczak et al. 2015). 
The brightest asteroids typically were targeted in one or more 
spectroscopic surveys— e.g., the Eight Color Asteroid Survey 
(ECAS, Zellner et al. 1985) or the Small Main-Belt Aster¬ 
oid Spectroscopic Surveys (SMASS; Xu et al. 1995, Bus & 
Binzel 2002). Dimmer objects however often only have color 
information from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000, Ivezic et al. 2001, Parker et al. 2008). The color 
index of Waszczak et al. (2015) puts asteroids of all sizes on a 
single, quantitative color scale (a proxy for spectral slope), the 
endmembers of which we identify with the C-type and S-type 
complexes. We use the terms ‘C types’ and ‘S types’ purely 
for compatibility with the literature, noting that our color in¬ 
dex combined with the classification thresholds we apply to it 
represent original definitions of these two groups. 

Use of a one-dimensional color metric sacrifices sensitiv¬ 
ity to secondary taxonomic groups such as M types, D types 
and V types. While our purpose is mainly to examine the 
first-order taxonomic variation, in certain cases we divide 
our color-defined ‘C types’ into two classes (Chigh and Q ow ) 
on the basis of visible and near-infrared albedo (as detailed 


in Section 3.2). At the end of this work we also examine 
well-established C-complex subgroups from both the Tholen 
and Bus/Binzel taxonomies (Section 6), showing how these 
known subgroups (e.g., M types and G types) are distributed 
in this work’s color-albedo space and highlighting these sub¬ 
groups’ unique UV properties. 

A second distinction between this work’s approach and 
Roettger and Buratti (1994) is that, rather than comparing the 
geometric albedo in the UV with that of the visible band, we 
focus on the difference in apparent magnitudes between the 
UV and visible. One motivation for doing this is we need not 
make any assumptions about the phase function of asteroids in 
the UV. The challenge however is that we must accurately es¬ 
timate the visible flux at the time of the UV observations. As 
discussed in Section 3, we adopt (and compare) two distinct 
methods for predicting the visual magnitude. The first method 
simply adopts the widely-used MPC 1 predicted magnitudes; 
the second method applies color-dependent phase-function 
and bond-albedo estimates adapted from the Waszczak et al. 
(2015) study of lightcurves from the Palomar Transient Fac¬ 
tory survey 2 (PTF; Law et al. 2009; Rau et al. 2009). 

2. GALEX ASTEROID OBSERVATIONS 

Extracting detections of known asteroids from a survey in¬ 
volves a three-dimensional (R.A., Dec., time) cross-matching 
of ephemerides against the survey’s time-stamped image 
boundaries (e.g., Ofek 2012). We modified software orig¬ 
inally used to search for asteroids in PTF (Waszczak et al. 
2013, 2015) to instead search for asteroids in GALEX. 

We first retrieved the metadata of all GALEX images, 
available from the Space Telescope Science Institute via 
command-line queries with the CasJobs tool (Li and Thakar 
2008). We then indexed all image centers with respect to 
(R.A, Dec.) into uniformly-spaced sky cells of 3-degree ra¬ 
dius. For all ~380,000 numbered asteroids, we queried JPL’s 
online service HORIZONS (Giorgini et al. 1996) to gener¬ 
ate a 1-day-spaced ephemeris spanning 2003-2012. Using an 
object-specific search radius equal to 3 degrees (cell radius) 
plus 0.75 degrees (FOV radius) plus the object’s maximum 
1-day motion (~ 10 arcminutes for most main-belt objects), 
we matched the ephemeris points against the sky cells. For 
each matched cell, we filtered out all images in that cell not 
within the epoch range of the matched ephemeris points, then 
for each surviving image we re-queried HORIZONS for the 
precise location at each observed epoch. We next performed 

1 IAU Minor Planet Center, http://minorplanetcenter.net 

- http://ptf.caltech.edu 
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Figure 2. Characteristics of positive asteroid detections from the two GALEX 

a 1.25-degree-radial match of these precise positions against 
the relevant GALEX image centers. 

We found ~850,000 predicted detections of numbered as¬ 
teroids (with no limit on apparent magnitude) in GALEX us¬ 
ing this method. For each predicted detection, using CasJobs 
we queried the GALEX single-visit source list (as opposed to 
the co-added source list). Multiple matches near the same 
point occurring more than 6 hours apart were excluded, as 
were all matches further than 2" from the predicted location. 
Additionally, to ensure the inclusion of greater than (approxi¬ 
mately) 5cr detections, we discarded all matches with NUV > 
21 mag in the shorter exposures (AIS program), and discarded 
all matches with NUV > 22.7 mag in the longer exposures 
(MIS program), following the limiting magnitudes quoted by 
Morissey et al. (2007). 

Following the above procedure and criteria, we extracted a 
total of 1,342 positive NUV detections of 405 unique aster¬ 
oids which were detected by GALEX at least twice (and no 
FUV detections, as expected). These detections are listed in 
Table 1; several histograms detailing these detections appear 
in Figure 2. 

All of the GALEX -observed asteroids are in the main-belt; 
the sample includes no near-Earth or outer-solar system ob¬ 
jects. 



position residual (arcsec) 



NUV mag 


surveys (distinguished by exposure time) shown separately in black and green. 

and precisely known at all observed epochs. Our two meth¬ 
ods for estimating V differ in their assumptions regarding 
(and observational data used to constrain) H and 0. In both 
cases we do not attempt to model the rotational term S , but 
rather incorporate <5 into the uncertainty of V using lightcurve 
amplitude estimates from the literature. In particular, 388 of 
the 405 GALEX -observed asteroids have an amplitude lower- 
limit estimate available in the Lightcurve Database (Warner et 
al. 2009, Harris et al. 2012). 

In the following sections we refer to two different albedo 
quantities. The visible-band geometric albedo py relates to 
the visible-band bond albedo Ab on d and the phase function 0 
(of Equation [1]) according to 

Pv = (/ 0(a) sin(a) da'j = (2) 

The above equation also defines the phase integral q. The 
bond albedo Ab on d is defined as the total visible light energy 
reflected or scattered by the asteroid (in all directions) divided 
by the total visible light energy incident upon the asteroid 
(from the Sun). Assuming the asteroid has a circular cross- 
section of diameter D, this can be expressed as 


3. MODELING VISIBLE MAGNITUDES 

In this section we consider two distinct methods of estimat¬ 
ing the visible magnitudes corresponding to all GALEX NUV 
detections; this in turn provides the distribution of the aster¬ 
oids’ NUV — V color. The general model for an asteroid’s 
apparent visual magnitude V (log flux) is 

V = H + (5 + 51og 10 (rA) - 2.5 log 10 [0(a)], (1) 

where H is the absolute magnitude (a constant), 5 is a periodic 
variability term due to rotation (e.g., if the object is spinning 
and has some asymmetry in shape or albedo), r and A are 
the heliocentric and geocentric distances (in AU), and 0 = 
0(a) is the phase function, which varies with the solar phase 
angle a (the Sun-asteroid-Earth angle). When a = 0 ( i.e. , at 
opposition), 0 = 1 by definition, while in general 0 < 0 < 1 
for a > 0 (with 0 decreasing as a increases). 

All asteroids for which we have extracted GALEX observa¬ 
tions have known orbits, meaning r, A, and a are accurately 


= £ /(a) sin (a)da 

b ° nd ~ (/ Su „/4^AU 2 ) x tt(D/2)2’ IJ! 

where /(a) = 10 - ' (“)/ 2 - 5 is the asteroid’s flux as a function 
of phase angle, with V (a) = H — 2.5 log 10 0(a) being Equa¬ 
tion (1) evaluated at S = 0 and r = A = 1 AU (similarly, 

/Sun = 10-W 2 -5). 

3.1. (H,G) from MPC data 

The first method for estimating V adopts the Minor Planet 
Center’s computed absolute magnitudes ( 11\m:A which are 
regularly updated by the MPC’s automated processes and uti¬ 
lize the Lumme-Bowell G-parameter model for 0 (Bowell et 
al. 1989). This same ( H , G) model then predicts the apparent 
magnitude I4ipc as a function of solar phase angle. 

The T/mpc values are fit to photometry provided by a va¬ 
riety of surveys/individuals, many of whom may use slightly 
different absolute calibration standards or filters with slightly 
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Figure 3. We compute Vptf model magnitudes by first assigning fixed Ab on d and G 12 values to each GALEX -observed asteroid depending on its color class; we 
then use D to compute H, and finally use the assumed G 12 value to predict V. The fixed values of Abond and G 12 are medians from the color-albedo-Gi 2 data 
in Waszczak et al. (2015), 2D histograms of which are shown here. Above each plot is the sample size (N = ...). A total of 642 asteroids have color data and 
G 12 values in the PTF data; 355 of these also have diameters available (required to compute Ab on d)- Panel D shows that WISE W 1 geometric albedos correlate 
with the PTF bond albedo; we thus use the WISE pwi data to assign C types either a low (Abond ~ 0.01) or high (Abond ~ 0.04) bond albedo. 


different specifications. A small fraction of asteroids have fit¬ 
ted G values; Harris & Young (1988) present mean G values 
for several major taxonomic classes, with G = 0.15 being 
an average between the C types (G ss 0.08) and the S types 
(G ~ 0.23). For the majority of asteroids the MPC uses an 
assumed G = 0.15 with this model. Waszczak et al. (2015) 
compares the E/mpc values with H magnitudes derived from 
a model that includes rotation and the more modern (II,Gi 2 ) 
phase function of Muinonen et al. (2010). Among bright as¬ 
teroids the relative difference is typically between 0.3% to 
3%, corresponding to (on average) an --0.07 mag discrepancy. 

Though H MP c values are available for all 405 GALEX- 
observed asteroids, we only consider the subset of 315 as¬ 
teroids having visible-band color indices of either less than 
0.25 (‘C types’) or greater than 0.75 (‘S types’). Of these, 41 
asteroids have Gmpc f =- 0.15. 

3 . 2 . ( D,Ab onl j,G 12 ) from PTF, infrared, and color data 

Our second means of estimating visual magnitudes applies 
only to asteroids having both a color index and a diameter 
estimate constrained from thermal fluxes in an infrared sur¬ 
vey 3 . In this approach we use the G 12 -parameter model for <j> 
(Muinonen et al. 2010), and we replace H with its equivalent 
expression 4 in terms of the diameter D, bond albedo Ah 
and phase integral q: 


*bond> 


H = —5 log 


O \/^'l !,on(i / q 


10 


(4) 


1329 km I 

where the phase integral q is a linear function of G 12 : 


( 0.2707 - 0.236G 12 ifG 12 <0.2; 
q (^i 2 ) j 0.2344 - 0.054Gi2 otherwise. 

We again define ‘C types’ as all asteroids with color indices 
less than 0.25 and ‘S types’ as all with color indices greater 
than 0.75. For S types we then consider diameters derived 
from any of four infrared surveys (see Footnote 3), while for 

3 Similar to the color data, the diameter data set we use is a compilation 
of products from several surveys and described in the appendix of Waszczak 
et al. (2015). The source IR surveys are WISE (Wright et al. 2010, Masiero 
et al. 2011, Masiero et al. 2014), IRAS (Matson et al. 1986, Tedesco et al. 
2002a), MSX (Tedesco et al. 2002b), and AKARI (Usui et al. 2011). 

4 Equation (4) follows directly from combining Equations (l)-(3). The 

constant 1329 km depends on somewhat arbitrarily-defined quantities such 
as the Sun's visual magnitude and the ratio of an AU to a kilometer. 


Table 2 

Abond an( -i C> 12 (based on PTF data) of color-defined taxonomic groups. 


class 

name 

color 

index 

WISE 

Pwi 

Abond 

median 

Abond 

scatter* 

G\2 

median 

G 12 

scatter* 

s 

> 0.75 

N.A. 

0.056 

0.016 

0.36 

0.16 

Chigh 

< 0.25 

> 0.125 

0.038 

0.022 

0.42 

0.20 

Clow 

< 0.25 

< 0.125 

0.010 

0.003 

0.84 

0.16 


‘Scatter is here defined as 0.5 X (84 th percentile —16 th percentile) 

C types we specifically require that the asteroid have been ob¬ 
served in the WISE 4-band cryogenic survey (Wright et al. 
2010, Masiero et al. 2014 and references therein). Both the 
WISE IV 1 -band geometric albedo pwi and the PTF-derived 
bond albedo 5 A holl( | show evidence of bimodality among ob¬ 
jects with color indices less than 0.25 (Figure 3 panel D). 
Thus, we divide the C types into low-bond-albedo (Ci ow ) and 
high-bond-albedo (Chi g h) subgroups based on their pwi as re¬ 
ported by Masiero et al. (2014). In Section 6 we show that 
the Chigh types most closely correspond to what other authors 
have called M types. 

Waszczak et al. (2015) computed Ab on d and G 12 values for 
— 1,600 asteroid lightcurves in the PTF survey. Using that 
work’s data (Figure 3) we compute median ,4 hon( | and G 12 
values (and associated scatter) for the S, Ci ow and Chigh tax¬ 
onomic groups. Table 2 summarizes the definitions and as¬ 
sumed Abond and G 12 values of these groups. There are 245 
GALEY-observed asteroids (out of the 405 in Table 1) which 
have color and diameter data available, allowing them to be 
modeled by this method. To each GALEY-observed asteroid 
we assign the appropriate Ab on d and G 12 value based on its 
class membership, then use its diameter to compute a model 
absolute magnitude (//pi h) using Equation (4). Together with 
the assumed G 12 value, this TTptf then predicts the apparent 
magnitude Vptf at each GALEY-observed solar phase angle. 

3 . 3 . Rotational uncertainty in V 

Both the U MPC and Vptf model magnitudes discussed here 
lack an estimate of the rotational term (5 in Equation [1]). We 
account for this by incorporating a term for rotational modu¬ 
lation into the reported uncertainty of V. Of the 315 asteroids 
with Vmpc values, 302 have an amplitude lower limit listed 
in the Lightcurve Database (Warner et al. 2009, Harris et al. 
2012), while for the 245 asteroids with Vptf predictions there 

3 The visible bond albedo Abond uses the same WISE diameter used by 
Masiero et al. (2014) in computing the W 1-band geometric albedo pwi ■ 
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Figure 4. Distribution of the NUV — V color for GALEX -observed asteroids using the H , G model with MPC data to predict V. Plots B and C define C types 
and S types as objects with color indices of <0.25 and >0.75, respectively. Plots C-F include only the subset with a diameter estimate available; this subset is 
precisely the same sample considered in Figure 6. 



Figure 5. Flowchart visualizing the steps in the GALEX -observed asteroid 
sample selection process. Each box is a subset of the box pointing to it. 


are 239 with reported amplitudes. As shown for instance by 
Waszczak et al. (2015), asteroids in the relevant size range 
typically have amplitudes less than ~0.4 mag. For the few 
objects in our sample lacking an amplitude limit, we assume 
a value of 0.2 mag. 

Assuming an asteroid’s rotational phase <p to be random at 
the time of a GALEX detection ( i.e ., with a probability dis¬ 


tribution of the form P« oc constant), then the probabil¬ 
ity distribution of a basic sinusoidal S (i.e., one of the form 
S = Sq sin ip) can be shown to have the form 


P(S) ex. 


1 


( 6 ) 


where So is the amplitude. We use Equation (6) as a proba¬ 
bility density function to generate, for each modeled V, a set 
of 10 4 simulated S values. These simulated S are added to an 
equal number of model V magnitudes computed by random 
(Gaussian distribution) sampling of the component terms: in 
the case of 14 ipc we just assume a fixed Hmpc uncertainty 
of 0.1 mag, whereas for the Vptf values we randomly sam¬ 
ple all three of Ab 0 nd> G 12 , and D, using the scatter values 
in Table 2 for the first two and the literature-reported diam¬ 
eter uncertainty for D. The 16 th to 84 th percentile spread in 
the distribution of combined S + V values then becomes the 
quoted uncertainty for V. 


4. NUV - V COLOR DISTRIBUTION 

Having computed the model V magnitudes, we obtain the 
NUV — V color for each GALEX asteroid detection and the 
corresponding uncertainty. The latter contains an additional 
rotational uncertainty component (now associated with the 
NUV observation), again determined by repeated sampling 
of Equation (4) as described above. Since all the asteroids 
we consider have more than one GALEX NUV detection, we 
compute the variance-weighted average NUV — V color for 
each asteroid (plotted in Figures 4 and 6); the uncertainty in 
this average is the inverse quadrature sum of the individual 
uncertainties. 

In Figures 4 and 6 (and the accompanying analysis) we have 
omitted all asteroids with NUV — V uncertainties of greater 
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Figure 6. Distribution of the NUV — V color for GALEX-observed asteroids using the //■ A|, nn( |.G'i 2 model with PTF data, infrared data, and color data to to 
predict V. See Table 2 for the definitions of the S, C] ow . and C^igh groups in plots B, C, E, and F. 


Table 3 

NUV — V color (mag. units) of GALES’-observed asteroids and sample sizes 


class 

name 

NUV 

median 

— Vmpc 
scatter* 

NUV 

median 

— Vptf 
scatter* 

Vmpc 

TVptf 

s 

6.52 

0.25 

6.71 

0.21 

72 

69 

c 

5.90 

0.19 

6.03 

0.22 

225 

154 

Chigh 

- 

- 

6.14 

0.33 

- 

29 

Clow 

- 

- 

6.02 

0.19 

- 

125 


‘Scatter is here defined as 0.5 X (84 th percentile —16 th percentile) 


than 0.5 mag. As a result, the sample size of asteroids with 
NUV— Vmpc estimates is 297 (out of the 315 quoted in Section 
3.1), while the sample with NUV — Vptf estimates is 223 (out 
of the 245 quoted in Section 3.2). Figure 5 graphically sum¬ 
marizes the sample selection criteria in a flowchart. In Fig¬ 
ures 4, 6 and 7, the errorbars on the color indices were com¬ 
puted by a bootstrapping process described in the appendix of 
Waszczak et al. (2015). 

Both the Vmpc and Vptf model magnitudes produce a bi- 
modal NUV — V color distribution, with the S types having 
the redder NUV — V color (panels A and B of both Figures 4 
and 6). Median and scatter of NUV — V for the various classes 
appear in Table 3. To formally ascertain the inequality of the 
two distributions, we use the two-sided Kolmogorov-Smimov 
(KS) test (Massey 1951), which compares two empirical dis¬ 
tributions via a bootstrap method. In particular this test com¬ 
putes a statistic quantifying the extent to which the cumula¬ 
tive distribution function differs in the two distributions being 
compared. For the Vmpc model we find the C-type NUV — V 
color distribution differs from that of the S-type distribution 
at an 11.6a significance level (Figure 4 panel B). For the Vptf 
model (Figure 6 panel B) we find the C types (Q ow and Chigh 
combined) differ from the S types at an 8. la level, while the 
Ciow and Chigh types only differ at a 1.9er level (this difference 


is thus not statistically significant). 

An important characteristic of our sample is that the C types 
outnumber the S types by a ratio of 3:1 in the Vmpc sample and 
a ratio of 2: 1 in the Vptf sample (cf. panel C of Figures 4 and 
6). This ratio is a combination of (1) the inherent difference in 
the population sizes of the two types (above a given diameter 
cut-off), a detection bias due to S types dominating the inner 
main-belt and thus typically having brighter apparent magni¬ 
tudes for a given size and albedo, and (3) the difference in the 
S and C types’ NUV albedo (discussed in Section 5). 

In panels C-F of Figure 4 the sample size decreases from 
N = 297 down to N = 223 asteroids as we consider only 
those objects in the Vmpc sample that also have available di¬ 
ameters (this is equivalently the Vptf sample considered in 
Figure 6). We compute the MPC-based visible bond albedo 
/1m pc using Equation (2) together with the asteroid’s TJmpc 
and Gmpc values. In particular, there are 38 asteroids (out of 
the 223 with diameters) with a measured Gmpc / 0.15; for 
the remainder we assume Gmpc = 0.15 for consistency with 
the manner in which the Uvipc are computed. Analogous to 
Equation (5), the phase integral for the G-model (required for 
computation of ,1 M pc via Equation [4]) is 

q(G) = 0.290 + 0.684G, (7) 

as given by Bowell et al. (1989). With the Gmpc = 0.15 
assumption for the majority of the asteroids in our sample, 
the /l mpc values are not expected to be as accurate as the 
,4bond values computed for instance by Waszczak et al. (2015), 
wherein distinct q values were fitted to each object on the ba¬ 
sis of a lightcurve. Nonetheless, it is instructive to compute 
41mpc, e.g., to check for consistency with the class-median 
Abond values, and to exploit as a second taxonomic metric in 
addition to visible color. 

Figure 4 panel E shows that NUV — V correlates with 

































































Asteroids in GALEX 


7 


N = 223 N = 223 N - 223 




Figure 7. Various checks for systematic differences in the predicted V magnitudes output by the two different photometric models. Bottom row: Investigation 
of phase-angle-dependence on the NUV — V color. 


A M pc (^Spearman = 0.698, >10cr significance), similar to 
how NUV — V correlates with the color index in panel A 
(^Spearman = 0.491, >10(7 significance). Unlike the color in¬ 
dex however, the separation between the Ci ow and C h i gh sub¬ 
groups is qualitatively evident in this plot. Figure 4 panel 
F combines all three parameters; note the axes are the same 
as Figure 3 panel B, with Ampc replacing Ab on d and the data 
consisting of GALEX- observed asteroids rather than PTF- 
observed asteroids. 

Figure 6 panel F confirms (independently of Figure 3 panel 
D) the validity of using WISE W 1-band geometric albedo as a 
proxy for visible bond albedo to separate Ci ow from Chigh—the 
two classes robustly differ in their Ampc distributions (9.5cr 
KS-test significance). However, the class-median Ampc val¬ 
ues of the Ciow, Chigh, and S types are 100%, 67%, and 63% 
greater than their class-median PTF-based /l hond values in Ta¬ 
ble 2. This reflects the differing values of H and q produced 
by the G and G 12 models, as well as the fact that we apply 
class-specific G12 values, whereas Gmpc = 0.15 is assumed 
for the majority of asteroids, regardless of their class. 

Consideration of both the 14ipc and Vptf model magnitudes 
provides two independent means of computing NUV — V; 
this helps rule out the effect of potential systematic errors 
unique to either one of the V models, as well as possible bi¬ 
ases in the distinct observational data sets upon which each 
V is based. In Figure 7 panels A-C we examine the distribu¬ 
tion of V M pc ^ Vptf for all 223 asteroids having both V esti¬ 
mates. The median of Vmpc — Vptf is 0.13 mag (scatter of 0.25 
mag), indicating the MPC-based model consistently produces 
brighter V estimates. For C types the median Vmpc ■ Vptf 
is 0.14 mag and for S types it is 0.12 mag; the two groups’ 
Vmpc — V>tf distributions differ with less than O.lcr signifi¬ 
cance in a KS-test. 

In Section 1 we motivated our choice to examine the dif¬ 
ference in apparent magnitude between UV and visible (as 


opposed to the difference in albedo in UV and visible) by not¬ 
ing that little is known of asteroid phase functions in the UV, 
rendering difficult the estimation of UV absolute magnitudes 
(and hence UV albedos). A potential issue with this approach 
which we have heretofore ignored is that, if the phase function 
does differ significantly in the UV from the visible, then the 
NUV — V color will vary with phase angle. Figure 7 panels 
D-F attempt to ascertain whether such a trend exists by con¬ 
sidering the median phase angle at which each asteroid was 
detected by GALEX. 

C types are observed at a median median phase angle of 7.2 
deg compared to the S types’ median median phase angle of 
9.0 deg. This is explained by the fact that C types on average 
have larger semi-major axes, which geometrically correspond 
to lower observed phase angles from Earth. Within the C-type 
group, median phase angle correlates with NUV — Vmpc at 
Pspearman = —0.1 (1.5(7 significance) and with NUV — Vptf at 
Pspearman = 0.01 (O.lcr significance). Among S types, median 
phase angle correlates with NUV — Vmpc at Pspearman = 0.01 
(O.lcr significance) and with NUV — Vptf at pspearman = 0.07 
(0.5(7 significance). We therefore cannot claim any phase an¬ 
gle dependence for NUV — V, regardless of the taxonomic 
group or V-model being considered. 

Various works (e.g., Sanchez et al. 2012 and references 
therein) discuss the phenomenon of asteroid phase reddening, 
i.e. an observed reddening of visible color with increasing 
phase angle. Very few survey-scale samples have been used 
to test for the presence of this effect. Szabo et al. (2007) com¬ 
puted slightly different phase-angle dependences for the g — r 
and r — i colors of Trojans in SDSS, though these relations 
were not separately computed for the Trojans’ two taxonomic 
groups. Waszczak et al. (2015) did not detect any statisti¬ 
cally significant difference between G 12 fits to r-band PTF 
lightcurves and g-band PTF lightcurves (among asteroids that 
had data in both bands). The extent to which a phase-function 
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dependence on wavelength exists between the UV and visible 
remains unclear. Future UV surveys such as ULTRASAT (Sa- 
giv et al. 2014) offer the most promising means of testing this 
hypothesis, especially because (unlike GALEX ) they will ob¬ 
tain sufficient numbers of observations to adequately sample 
UV lightcurves, thereby providing the best possible data set 
for fitting UV phase functions. 

5. ALBEDO VS. WAVELENGTH 

If we assume that the phase function does not differ sig¬ 
nificantly between the UV and visible (or take Figure 7 pan¬ 
els E and F as justification of this statement), then we can 
compare the relative bond albedo versus wavelength for the 
different taxonomic groups using measured colors, filter re¬ 
sponse functions and the solar spectrum. If the phase func¬ 
tion does in fact vary significantly with wavelength, then this 
approach only provides the relative geometric albedo versus 
wavelength (see Equations [2] and [3]). 

Assume photometry from two filters (1 and 2) produce the 
color measurement ni \ — m 2 . This color relates to the solar 
flux distribution S( A), the albedos in each band (A \ and A 2 ) 
and the filter responses .Fj(A) and F 2 ( A) according to 

m(mi - m2 )/2. 5 = fF 1 (X)X~ 2 dXfS(X)F 2 (X)A 2 dX 

f S(X)F 1 (X)A 1 dX f F 2 (X)X~ 2 dX ( ’ 

which we adapted from a similar equation in Pickles (1998). 

Using NUV as band 1 and V as band 2 in Equation 6, we 
use the colors in Table 3 (specifically, the VpxF-based col¬ 
ors) to obtain the albedo ratio ,4 NUV /Ay, with uncertainties 
coming from the associated scatter in the colors. In Figure 
7 we plot these albedo ratios for the C types and S types, 
incorporating an additional uncertainty component from the 
transformation from r to V (see Waszczak et al. 2015 for a 
discussion of this transformation in the context of asteroids). 
The end-computed values are (A.Nuv/v4 r )c = 0.63and 
(A.Nuv,Mr)s = 0.33j)g Qg. The relative albedo values in the 
SDSS bands included for comparison in Figure 8 are taken di¬ 
rectly from a figure in Ivezic et al. (2001); likewise the ECAS 
data are taken directly from a figure in Zellner et al. (1985). 

Note that Figures 8 and 9 ignore the fact that the S types’ 
and C types’ absolute albedo in r-band differs. In other 
words, these plots could be converted into ones with abso¬ 
lute albedo on the vertical scale by multiplying the blue and 
red lines by their respective absolute r-band albedos, which 
would be similar to those listed in Table 2 for V band. As¬ 
suming (from Table 2) that C types have py = 0.054 ±0.019, 
then (pnuv)c = 0.035 ± 0.019, and assuming S types have 
py = 0.264±0.031, then (Anuv)s = 0.088 ±0.028. That is, 
S types have a higher mean NUV albedo than C types, similar 
to the visible. 

Both C and S types show a continued trend of decreasing 
albedo at shorter wavelengths. Whereas for S types this be¬ 
havior was already well-established in the 300-800 nm re¬ 
gion, for C types the u — g color had previously represented a 
significant deviation from the shallower slope observed from 
400-800 nm. The C-type NUV albedo in Figure 8 confirms 
the presence of a marked drop in albedo somewhere in the 
200-400 nm range. Given the resolution of Figure 8 and 
the uncertainties in the data points, we cannot judge whether 
the C-type albedo levels off between NUV and u band, or 
whether the slope between u and g bands persists into these 
shorter wavelengths. The IUE spectra from Roettger and Bu- 
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Figure 8. Top : Relative albedo (each type’s r-band albedo normalized to 
unity) versus wavelength for C types and S types. The leftmost (NUV) points 
are computed from this work’s data, the remaining albedos (SDSS bands) are 
taken directly from Ivezic et al. (2001). Bottom'. Bandpass response functions 
(colored lines) multiplied by the solar spectrum (black line). 



Figure 9. Top'. GALEX, SDSS, and ECAS C-type data from Figure 8 com¬ 
pared to //.ST-derived albedos for 1 Ceres and 21 Lutetia. Note the red color 
here has a different meaning than it does in Figure 8. 

ratti (1994) indicate C-type albedo is constant at least in the 
range 240-300 nm (as does the Lutetia data described below), 
so that the former interpretation may be more accurate. 

6. COMPARISON TO 11ST DATA 
6.1. Lutetia 

Weaver et al. (2010) obtained HST photometry of asteroid 
21 Lutetia in UV and visible bands. Lutetia has been classi¬ 
fied by various authors as an M-type asteroid; in the context 
of this work its color index is 0.05 (making it a C type) and its 
AImpc = 0.06 suggest it to be a Chi g h type in particular, though 
in this work’s system we formally would require a WISE pwt 
measurement to classify it as such. In the following section 
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we show that M types (a group in the Tholen taxonomic sys¬ 
tem) and our Chigh types are largely the same population. 

The HST Lutetia photometry revealed a steep drop in 
albedo around ~300 nm and nearly constant albedo in the 
200-300 nm region at a factor ~0.6 times the visible (r-band- 
equivalent) albedo. The HST observations of Lutetia thus 
generally agree with the C-type albedo trend (Figure 9), the 
main difference being the location of the UV albedo drop-off 
(the bluest two ECAS bands also demonstrate this difference 
between M types and C types, e.g., see Figure 2 of Bus et al. 
2002). The Rosetta spacecraft’s flyby of Lutetia enabled FUV 
observations with the on-board Alice UV imaging spectro¬ 
graph (Stern et al. 2011); the longest wavelengths of the FUV 
data (~ 190 nm) yield an albedo consistent with the constant 
value measured in the 200-300 nm range by HST, namely 
~0.6 times that of the visible albedo. 

6.2. Ceres 

HST photometry of asteroid 1 Ceres has also been obtained 
in the UV and visible (Parker et al. 2002, Li et al. 2006). With 
a color index of 0.01, Ceres is also a C type in our classifica¬ 
tion scheme, though its Ampc = 0.033 makes its placement in 
our Ciow vs. Chigh groups ambiguous (see Figure 6 panel F). 
Like Lutetia, Ceres lacks a reported pw i so that we cannot 
formally classify it as either Ci ow or Chigh- 

Ceres was observed by GALEX and thus is included in 
our MPC-data-based analysis; our measured NUV — Vmpc = 
6.45 ± 0.19 for Ceres make it a clear outlier from the C-type 
NUV — Vmpc distribution (Figure 4 panel B). In the Tholen 
taxonomic system Ceres is classified as a G type; in the fol¬ 
lowing section we show that other G types exhibit similarly 
high NUV — I 'mpc; values but less anomalous NUV — Vptf- 

The Parker et al. (2002) HST data show that around ~300 
nm Ceres’ albedo drops to as low as ~0.3 times the visible- 
band albedo—compared to the factor of ~0.6 seen for GALEX 
C types and the Lutetia data—but that around ~200 nm it 
appears to rise again to a more typical C-type UV albedo. 
Roettger and Buratti (1994) did not observe this unusually 
deep absorption feature near 300 nm in their IUE spectrum of 
Ceres; if real this feature could partially explain the anoma¬ 
lous NUV — Vmpc we observe for G types in GALEX. Figure 
9 shows Ceres data in the three HST bands observed by Li et 
al. (2006), none of which sample the 300-nm region contain¬ 
ing the putative absorption band, though these three bands do 
generally match the GALEX C-type data. 

7. C-TYPE SUBGROUPS 

C types deserve further consideration for several reasons: 
(1) C types outnumber S types in the GALEX samples by a 
factor of several, (2) our division of C types into Ci ow and 
Chigh merits interpretation in more conventional taxonomic 
systems, and (3) both of the //ST-observed asteroids in the 
previous section are known members of C-type subgroups, 
the UV properties of which are worth confirming with addi¬ 
tional group members. 

Figures 10 and 11 detail the distribution of GALEX- 
observed asteroids belonging to six classes each from the 
Tholen and Bus/Binzel taxonomic systems (Tholen 1989; Bus 
& Binzel 2002), the latter is sometimes referred to as the 
SMASSII system after the survey data with which it was de¬ 
rived. These two classification systems were created on the 
basis of different visible-band color data; a comparison of 
their group definitions is given in Table 1 of Bus et al. (2002). 
We consider only the subset of GA LEX-observed asteroids 


having both Vmpc and Vptf model magnitudes and omit sub¬ 
groups containing less than three objects. In the following 
subsections we briefly comment on these subgroups. 

One key interpretation of these data—supported also by the 
HST data in Figure 9—is that NUV-band albedo is not very 
useful for discriminating C-type subgroups, e.g., M types ver¬ 
sus G types, whereas u band appears to be more diagnostic 
in this regard. The u-bantl discrepancy between these sub¬ 
groups was remarked most notably by Zellner et al. (1985) 
in the ECAS data, but it was unknown at that time (indeed, 
up until now) whether the discrepancy in UV albedo became 
more or less pronounced shortward of ~300 nm. The NUV 
data indicate that the discrepancy lessens in the NUV, as M 
types do in fact exhibit a step-down in albedo (between NUV 
and u bands), similar to the step down the G types exhibit 
within u band. 

7.1. X complex 

The Tholen system’s X-type group includes asteroids with 
relatively flat visible color, including no substantial absorp¬ 
tion in the blue (in contrast to, e.g., the u-band drop-off seen 
in G types). The subgroups within the X group include E, M, 
and P types and are distinguishable only by albedo. 

The twelve M types in our sample all have Ampc > 0.03 
and pw i > 0.125, the latter formally makes them all Chigh 
types in this work’s classification system. The M types have 

NUV — Vmpc = 5.89 ± 0.15 and NUV - U PT f = 6.12 ±0.21, 

neither of which significantly differ from the C-type averages 
given in Table 3. This is consistent with the above-noted ob¬ 
servation that M-type Lutetia’s NUV — V is similar to that 
of the GALEX C types, despite an obvious difference in ri¬ 
band (Figure 9). Assuming all 29 of the Chigh types in the 
GALEX sample are in fact M types, then the Chigh types’ 
slightly higher NUV — Vptf = 6.14 ± 0.33 (compared to 
NUV — Vptf = 6.03 ± 0.22 for the whole C type group) 
agrees well with the M types’ slightly higher average. 

Complementary to the M types, the eleven P types in our 
sample all have Ampc < 0.03 and pw t < 0.125, the latter 
formally makes them all Q ow types. The P types have NUV — 
Vmpc = 5.74 ± 0.17 and NUV - U PTF = 5.91 ± 0.13. These 
values are less than both models’ C-type averages as well as 
less than the Q ow average, suggesting our Q ow group includes 
more diverse objects than just P types (e.g., the five F types 
also all have Ampc consistent with Ci ow ). 

There are 14 GALEX- observed asteroids listed simply as X 
types in the Tholen system (presumably because no visible 
albedo was available at the time of classification); Figure 10 
shows that these are in fact distributed across both the Q ow 
and Chigh albedo ranges. 

In the Bus/Binzel system, the X complex consists of four 
subgroups: Xc, Xk, X and Xe, these being differentiated by 
their spectral slope and presence of various absorption fea¬ 
tures. In the GALEX sample the most numerous of these are 
the Xc types, which have the least red visible color and seem 
to include both high and low visible albedo members. Both 
the Xe and Xk types have higher visible color indices (with 
larger uncertainties in the color). As with the Tholen X types, 
we see no systematic trends with respect to the NUV proper¬ 
ties of these subgroups. 

7.2. G types 

Three GALEX- observed asteroids are categorized as G 
types. Like G-type Ceres, these have intermediate Ampc and 
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Figure 10. Visible color/albedo distributions of Tholen-classified (Tholen 1989) and Bus/Binzel-classified (Bus & Binzel 2002) C-type subgroups among this 
work’s sample of GALEX -observed C-type asteroids. See text for further information. 
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Figure 11. NUV — V color distributions of Tholen-classified (Tholen 1989) and Bus/Binzel-classified (Bus & Binzel 2002) C-type subgroups among this work’s 
sample of GALEX -observed C-type asteroids. See text for further information. 


































12 


Waszczak et al. 


an above-average NUV — V'vipc = 6.22 ± 0.11. In contrast, 
however, the G-type NUV — Vptf = 5.79 ± 0.10 lies slightly 
below the C-type average. The reason for this discrepancy 
is that all three G types in this sample have pw\ < 0.125 
and so are formally classed as Ci ow objects, as a result their 
assumed Hb on d = 0-01 in the computation of Vptf may be 
too low. On the other hand, Oszkiewicz et al. (2011) fit 
G 12 = 0.88 ± 0.2 to Ceres’ phase function, suggesting that 
the assumed G 12 = 0.84 ± 0.14 for Ci ow types (Table 2) is 
a more valid assumption for G types than the Chigh value of 
G 12 = 0.42 ± 0.20. Hence the G types seem not to fit well 
into either of our Ci ow or Chi g h groups, and hence are not ac¬ 
curately modeled by our Vptf- 

The three-asteroid G-type sample’s higher than average 
NUV — Vmpc agrees with the Ceres HST data (Figure 9), 
which as discussed above could be indicative of an absorp¬ 
tion feature at ~300 nm unique to G types (Li et al. 2006), 
the precise shape and location of which remains unresolved 
in the broadband photometry considered here. 

Tholen’s G types are represented in the Bus/Binzel sys¬ 
tem by the Cg and Cgh groups; however no asteroids in our 
GALEX sample have either of these SMASSII labels. 

7.3. B types 

Members of the Tholen B and F classes, represented also 
by the Bus/Binzel B and Cb classes, all are classified as Ci ow 
types in the GALEX sample based on their pwi- Unlike the G 
types, the B types are not anomalous in NUV — U M pc, mean¬ 
ing the B types likely lack the G types’ strong absorption at 
300nm. The B types also are characterized by slightly higher 
Hmpc = 0.026 compared to the Ci ow average Ampc = 0.020. 
Hence, like the G types, the B types show a lower than aver¬ 
age NUV — Vptf symptomatic of an underestimated /l hon d and 
therefore too dim of a predicted Vptf. 

8. SUMMARY 

We present NUV-band photometry of 405 asteroids ob¬ 
served serendipitously by GALEX from 2003-2012. Using a 
compilation of visible-band color data, we select the subset of 
these GALEV-observed asteroids belonging to the C-type or 
S-type classes. We then compute the visual-band magnitude 
(using two different models) corresponding to each GALEX 
detection in an effort to study the NUV — V color. For both V 
models, the derived NUV — V color distribution is bimodal, 
with S types having the redder color, just as they do within 
the visible band. The average C-type NUV — V agrees with 
HST observations of the asteroids Lutetia and Ceres, both of 
which are members of the visible-color-defined C-type group. 
Slight differences in the measured NUV — V among known 
taxonomic subgroups of the C types may indicate membership 
in either the M-type or G-type subgroups, though the 300-400 
nm region ((/-band) is more diagnostic of this division. 
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GLOSSARY OF ACRONYMS AND SYMBOLS 

For the reader’s convenience. Table 4 summarizes the vari¬ 
ous acronyms and mathematical symbols used in this work. 
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Table 4 

Glossary of acronyms and symbols used repeatedly in this work. 


symbol or 
acronym 
a 

Abond 

Ampc 
Anuv 
Chigh types 
Ciow types 
D 

ECAS 

FUV 

G 

G12 
G types 
GALEX 
H 
HST 
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V 
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Minor Planet Center, http://minorplanetcenter.net 
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